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SUMMARY 


An effective computational strategy is presented for the large-rotation, non- 
linear, axisymmetric analysis of shells of revolution- A total Lagrangian descrip- 
tion of the shell deformation is used, and the analytical formulation is based on a 
form of Reissner's large-deformation theory including the effects of transverse shear 
deformation, laminated orthotropic material response, and moments acting about the 
normal to the middle surface. Only axisymmetric deformations are considered, and the 
fundamental unknowns consist of six stress resultants and three generalized displace- 
ments of the shell. The element characteristic arrays are obtained by using the 
Hellinger-Reissner mixed variational principle. The polynomial interpolation (or 
shape) functions used for approximating the stress resultants are of lower degree 
than those used for approximating the generalized displacements. 

The three key components of the proposed computational strategy are as follows: 

(1) use of mixed finite-element models with discontinuous stress resultants at the 
element interfaces, thereby allowing the elimination of these stress resultants on 
the element level; (2) substantial reduction in the total number of degrees of freedom 
through the use of a multiple-parameter reduction technique based on perturbation 
expansion; and (3) reduction in the size of the analysis model through the decomposi- 
tion of asymmetric loads into symmetric and antisymmetric components coupled with the 
use of the multiple-parameter reduction technique. 

The potential of the proposed computational strategy is discussed. Numerical 
results are presented to demonstrate the high accuracy of the mixed models developed 
and to show the effectiveness of using the proposed computational strategy with these 
models . 


INTRODUCTION 

In recent years nonlinear analysis of static and dynamic problems has become the 
focus of intense research efforts. The increasing importance of nonlinear analysis 
is largely due to the emphasis placed by manufacturers, contractors, and certifying 
agencies on realistic modeling and accurate analysis of critical structural compo- 
nents. This analysis endeavor has prompted the development of versatile and powerful 
finite-element discretization methods as well as development of improved numerical 
methods and programming techniques for nonlinear analysis of structures. As this 
nonlinear analysis technology becomes more sophisticated, it can be used for pre- 
dicting the response of a wider class of structures. One of the challenging non- 
linear problems is that of predicting the response of aircraft tires during ground 
handling operations. The tires must be able to taxi for long distances under heavy 
loadings, to operate at speeds greater than 200 knots, to absorb much of the energy 
of the landing impacts, and to provide steering and braking forces for safe ground 
handling operations. Yet these tires must be as lightweight and durable as possible. 
To date, only a few analytical tools are available to assist the tire designer, and 
current design technology of tires is based largely on empirical methods. A summary 
of past and current modeling techniques for tires is presented in reference 1. 

In addition to the diverse and severe loading conditions to which aircraft tires 
are subjected, several other factors contribute to the difficulty of constructing 



analytical models of tires. The tire is a composite structure composed of rubber- 
textile constituents which exhibit anisotropic, nonhomogeneous material properties. 
The ground handling operations routinely result in very large rotations and large 
deformations in the tire. Also, the tire carcass is thick enough to allow signifi- 
cant transverse shear deformation to occur. All these attributes work, in concert to 
make the development of practical analytical tire design tools extremely difficult. 

A shell finite-element model incorporating the effects of large rotations, laminated 
anisotropic material response, and transverse shear deformations, however, shows 
promise for tire modeling applications. 

The present study constitutes a first step towards the development of computa- 
tional models for tires. The objectives of this paper are twofold: (1) to present 

an accurate and efficient computational strategy for the nonlinear, axisymmetric 
analysis of shells of revolution, and (2) to demonstrate the effectiveness of the 
proposed strategy by means of numerical examples. The three key components of the 
proposed computational strategy are as follows: (1) use of mixed finite-element 

models with discontinuous stress resultants at the element interfaces; (2) substan- 
tial reduction in the total number of degrees of freedom through the use of a 
multiple-parameter reduction technique based on perturbation expansion; and 
(3) reduction in the size of the analysis model for symmetric structures through the 
decomposition of asymmetric loadings into symmetric and antisymmetric components 
coupled with the use of the multiple-parameter reduction technique. To sharpen the 
focus of the present study, axisymmetric static and large-rotation problems of shells 
of revolution are considered. A total Lagrangian description of the shell deforma- 
tion is used, and the analytical formulation is based on a form of Reissner’s large- 
deformation theory with the effects of transverse shear deformation, laminated 
orthotropic material response, and moments acting around the normal to the middle 
surface included. The mixed models presented herein have both stress resultants and 
generalized displacements as fundamental unknowns. The stress resultants are allowed 
to be discontinuous at interelement boundaries. 


MATHEMATICAL FORMULATION 
Variational Principle 

The analytical formulation is based on a form of Reissner's large-deformation 
theory of shells of revolution (see ref. 2 and appendix A) with the effects of trans- 
verse shear deformation, laminated orthotropic material response, and moments acting 
around the normal to the middle surface included. Only axisymmetric deformations of 
the shell are considered, and the load can be displacement dependent as well as non- 
conservative. The load is represented by two independent parameters q(^) and q(2) 
A mixed formulation is used in which the fundamental unknowns consist of the stress 
resultants N s , Nq , M s , Mq , Q s , and M n and the generalized displacements of the 
middle surface u, w, and <J> . (See fig. 1 for sign convention.) A total Lagrangian 

description of the shell deformation is used, and the shell configurations at differ- 
ent load levels are referred to the initial coordinate system of the undeformed 
shell . 

The mixed variational principle used in the element development is given by 


6n(N s ,N 0 ,M s ,M 0 ,Q s ,M n ,u,w,({)) 


6K - 6W = 0 


( 1 ) 
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where 
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In equations (2) and (3) , the quantity <J> 0 is the angle between the axis of revolu- 
tion and a normal to the undeformed middle surface; <J> is the rotation of the shell 
middle surface; r is the radius of the undeformed parallel circle; a, b, and g 
are shell compliance coefficients (inverses of shell stiffnesses; see appendix A); 

P u and P w are the intensities of external distributed loads in the radial and 
axial directions; m is the intensity of external distributed meridional moment; 
v P w' and M are the external concentrated loads and meridional moment (or total 

values of ring loads and moments); and 9 = — . The quantities in equation (3) with 

a tilde (~) denote prescribed forces and displacements. The summation sign in equa- 
tion (3) extends over all applied external forces and prescribed displacements in the 
shell meridian. 


Configuration-Dependent Load 


For the case of uniform hydrostatic pressure, the direction of the load always 
remains perpendicular to the deformed middle surface of the shell; thus, the load is 


3 



configuration dependent. In this case, the expressions for the external load 
components are given by 


P u = P o (1 + e so } i 1 + 7) Bin{ * + " V 

P w = P o (1 + £ so> i 1 + 7) cos <* + ' V 

m = 0 


( 4 ) 

( 5 ) 

(6) 


where p Q is the pressure load (per unit area of the deformed middle surface) ; e 
is the extensional strain in the meridional direction; and y Q is the transverse 
shearing strain. (See appendix A.) 


Finite-Element Discretization 

The finite-element discretization is performed by dividing the shell meridian 
into finite elements. The generalized displacements u, w, and <(> are approximated 
by sums of products of shape functions W and of nodal displacements X i . The 
stress resultants N g , Ng , M g , Mg, Q g , and M n are approximated by sums of pro- 
ducts of shape functions R and stress-resultant parameters H_. 

J 

The governing finite-element equations for each individual . element are obtained 
by first replacing the generalized displacements and stress resultants with their 
expressions in terms of the shape functions and then applying the variational 
principle (eqs. (1) to (3)). If the nodal displacements Xj and stress-resultant 
parameters Hj are varied independently and simultaneously, one obtains the follow- 
ing set of equations for each element: 
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is defined by the following equation: 


where the array 





(a) 



x. 


(9) 


The F and S terms are the linear flexibility and strain-displacement coefficients; 

and are normalized consistent-load coefficients; and are 

normalized load stiffness coefficients, which are symmetric for conservative loadings 

and asymmetric for nonconservative loadings; q^ and q^ are load parameters; 

A, B, A, and B are multipliers of the trigonometric contributions; Rj are 

integrals of shape functions (which cancel with the Aj°^ terms when <J> ^ = o); and 
superscript a denotes that the quantity is evaluated at the numerical quadrature 
point. The summation sign extends over the n numerical quadrature points in the 
element. The range of the lowercase indices is 1 to 3m (where m is the number of 
displacement nodes) , and the range of the uppercase indices is 1 to 64 (where 4 
is the number of nodes or parameters used in approximating each of the stress 
resultants within the individual elements) ; a repeated subscript in the same term 
denotes summation over its full range. The explicit forms of the arrays F, S, 

A, B, A, B, R, T, P^, P^ 2 \ Q^, and are given in appendix B. 

Note that in equation (8) the two loads P^ + Qj_j x j are assumed to be linear 
functions of the nodal displacements Xj (which is a reasonable approximation for 
the case of uniform hydrostatic pressure in eqs . (4) to (6)). The validity of this 
approximation is discussed in the section entitled "Numerical Results." Other 
displacement-dependent and nonconservative loads are discussed in reference 3. 

The continuity of the stress resultants N s , , M g , Mq , Q s , and M R is 
not imposed at interelement boundaries. Therefore, the stress-resultant param- 
eters Hj can be eliminated on the element level from equations (7) and (8) , and 
one obtains the following elemental equations in the nodal displacements X . : 
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The arrays A, B, C, D, £, and 7 are obtained from the definitions of the 
sponding_arrays A, B, C, D, E, and 7 by replacing the respective arrays 
B, A, B, and T in equations (12) to (21) with A, B, A, B, and -T. Sub 


(11) 

( 12 ) 

(13) 

(14) 

(15) 
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expressions which occur repeatedly in equations (11) to (21) are underlined with 
dashed, dotted, or dashed -dotted lines as appropriate. The linear finite-element 
equations for an individual element can be obtained from equation (10) by setting 


sin 





X . 
D 


COS <J) 


(a) 


1 


and neglecting the nonlinear terms. The resulting equations can be cast in the 
following form: 


(k. . - 

\ ID 


q (1) Q- (1) -q (2) Q 


( 2 )\ 

i-j / 


X. - q (1) pf X) + q (2) P. (2) 
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( 22 ) 


* 

where are the linear stiffness coefficients (excluding the load stiffness con- 

tributions Qij) for individual elements given by 


(e. (a) + A. (a) T + E E (£. (a ' 3) + D. (a ' e) Tf 6) 
\ ID 1 : / a e \ ^ 1 D 


(23) 


The terms are absent in the linear equation (22) because they cancel out the 

combination of Z and Z £ 

a 1 a 3 1 


Equivalence Between Mixed Models and Displacement Models 

A fixed model with discontinuous stress resultants at interelement boundaries is 
said to be zquuvcULdVVt to a displacement model if, after elimination of the stress- 
resultant parameters (on the element level) , the resulting finite-element equations 
for the two models are i. dentlcdt (i.e., eq. (10) for the mixed model is identical to 
the corresponding equation for the displacement model) . The models are said to be 
ndOJlty zqLLtVCLZ&ftt if their finite-element equations are almost LdzwtidCit (e.g., if 
the corresponding numerical coefficients in both sets of equations differ by a small 
percentage) . 

The equivalence between mixed models and some classes of displacement models was 
first established in reference 4 for linear problems. It was extended to moderate- 
rotation nonlinear problems of curved beams and shallow shells in references 5 and 6. 

For axisymmetric , large-rotation, nonlinear problems of shells of revolution 
with constant material properties, the following two groups of equivalent mixed and 
displacement models can be identified: 

Group 1 - Equivalent mixed models and standard displacement models . The condi- 
tions for the equivalence are as follows: 


7 



!• The same number of nodes m and Lagrangian interpolation functions are 
used for approximating the generalized displacements u, w, and <j> in 
both the displacement and mixed models. 

2. The number of quadrature points n in the Gauss-Legendre quadrature 
formula used for the numerical evaluation of all terms of the mixed and 
displacement models is equal to the number of displacement nodes. 

3. The polynomial interpolation functions for the stress resultants N s , 

V M s , Mq, Q s , and M n in the mixed models are of the same degree as 
those of displacements and the stress resultants are allowed to be dis- 
continuous at interelement boundaries. 

Group 2 - Equivalent mixed models and reduced integration displacement models . 

The conditions for the equivalence are as follows; 

1 • The same number of nodes m and Lagrangian interpolation functions are 
used for approximating the generalized displacements u, w, and (J> in 
both the displacement and mixed models. 

2. The same low-order Gauss -Legendre quadrature formula (i.e., n < m) is 
used for the numerical evaluation of all integrals in the displacement 
and mixed models. 

3. The interpolation functions for the stress resultants N g , N 0 , M g , M 0 , 

Q s , and M n in the mixed models are polynomials of lower degree than 
those of the displacement models, and the stress resultants are allowed 
to be discontinuous at interelement boundaries. The number of stress- 
resultant interpolation functions 4 equals the number of quadrature 
points per element n. Of particular interest is the case where 

n = 4 = m - 1 (where m is the number of displacement nodes) , and thus 
the interpolation functions for the stress resultants are polynomials 
one degree lower than those of displacements. 

Henceforth, the mixed and displacement models of group 1 are designated by MDm-m 
and DEm. The corresponding models of group 2 are designated by MDm-n and DRm-n. The 
aforementioned equivalences are based on comparing the dXdct dnotytld dXyKZkbAJOYl 6 ^OK. 
£,hd of the mixed and displacement models. Numerical experi- 

ments have demonstrated the superior performance of the models of group 2 over those 
of group 1. This result is demonstrated in the "Numerical Results" section. 


Solution of Finite-Element Equations 


The assembled finite-element equations (obtained from eqs. 
shell can be represented in the following compact form: 


(7) and (8)) for the 
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(24) 
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where [F] is a block diagonal matrix containing the elemental contributions f jj/ 
[$] is a rectangular matrix containing the Sjj contributions; {©(X)} and 

{flt(H,X)} are the vectors of trigonometric contributions; {x} is the assembled 
vector of nodal displacements; {H} is the assembled vector of stress-resultant 
parameters; {R} = {(5(0)} is a constant vector composed of the elemental contribu- 
tions Rj; (X) } and (X)} are the vectors of normalized external forces; 

and superscript t denotes transposition. Note that the linear global generalized 
stiffness matrix is given by 


-F 



Symm. 0 


Symm. 








=0 


(25) 


where I ranges from 1 to 6-6 (total number of stress-resultant parameters in the 
discretized shell) ; and i and j range from 1 to 3m (total number of displace- 
ment nodes in the discretized shell) . 


The solution of equation (24) is usually carried out by using an incremental- 
iterative technique (i.e., a predictor-corrector continuation method) wherein the 
vectors {h} and {x} corresponding to a particular pair of values of the load 

parameters q^ and q^ are used to calculate suitable approximations 

(predictors) for {H} and {x} at different values of q^ and q^. These 
approximations are then chosen as initial estimates for { H } and {x} in a 
corrective-iterative scheme such as the Newton-Raphson technique. 


MULTIPLE-PARAMETER REDUCTION TECHNIQUES 
Degree s-of -Freedom Reduction 

In recent years considerable progress has been made in the development of tech- 
niques for reducing the total number of degrees of freedom of the discretized struc- 
ture without sacrificing the accuracy of predicting the nonlinear response. Such 
techniques are called reduction methods and are essentially hybrid procedures com- 
bining the finite-element method with the classic Rayleigh-Ritz technique. (See 
refs. 7, 8, and 9.) The effectiveness of reduction methods for the nonlinear and 
postbuckling analyses of structures subjected to multiple independent loads has been 
demonstrated in references 10 and 11. A possible approach for applying multiple- 
parameter reduction methods in conjunction with the mixed formulation for the large- 
rotation, nonlinear analysis of shells of revolution is to use finite elements for 
the initial discretization of the shell and to express the vectors of unknown stress 
resultants {H} and generalized displacements {X} as a linear combination of a 
small number of global approximation (or basis) vectors. A Rayleigh-Ritz (or Bubnov- 
Galerkin) technique is then used to approximate the finite-element equations by a 
reduced system of nonlinear equations. This approach was applied to nonlinear large- 
rotation problems of curved beams in reference 12. 
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An effective choice for the basis vectors consists of a solution of the nonlinear 
finite-element equations (eqs. (24)) and its various-order path derivatives. Both the 
solution and the path derivatives are obtained at a particular pair of values q(D 

and q^ 2 t For the case of two independent load parameters q^ and q^, the 
path derivatives are given by 


{H} = [r] {^} 
{X} = [r] {<(,} 


where 


[r] = 



(26) 


(27) 


[r] - 



(28) 


In equations (27) and (28) , A^ an< ^ X 2 are parameters and {\j>} is the 

vector of undetermined coefficients, which are functions of X l and X 2 - The number 
of basis vectors H. is considerably smaller than the total number of degrees of 
freedom of the shell model. The path parameters Aj and A 2 can be identified with 
load, displacement, or pseudo arc-length parameters. (See ref. 11.) 


The equations used in evaluating the basis vectors are obtained by successive 
differentiation (with respect to A^ and A 2 ) of the finite-element equations of the 
discretized shell (eqs. (24)). The explicit forms of the first few of these equa- 
tions are given in appendix C. The path derivatives of the stress resultants can be 
eliminated on the element level, thereby reducing the size of the equations required 
for evaluating the displacement path derivatives. Note that even for the noncon- 
servative load case, the equations used in the generation of the initial basis 
vectors are symmetric. (See appendix C.) 


The explicit forms of the reduced system of equations in terms of the unknowns 
{i(j} are given in appendix D. Note that for large-rotation problems, the form of 
these equations is considerably more complicated than the corresponding equations for 
moderate rotations. In both cases, however, the number of reduced degrees of freedom 
is small. The computational procedure used in conjunction with the reduction method 
is described in references 11 and 12 and is not repeated herein. 


Model-Size Reduction Through Symmetry 

It has long been recognized that the size of the analysis model of a structure 
can be reduced by exploiting the d pfLLOSvi known symmetries of the response of that 
structure. Considerable effort has been devoted to identifying symmetries exhibited 
by various types of structures and exploiting these symmetries in their finite-element 
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analysis. For linear symmetric structures subjected to asymmetric external load, 
symmetry concepts can be used by decomposing the load into symmetric and antisymmetric 
systems. It is generally assumed, however, that such a decomposition is not useful 
for nonlinear problems in which the principle of superposition is not applicable. In 
fact, the nonlinear response to an antisymmetric load is, in general, asymmetric. In 
this section a simple procedure is outlined for reducing the size of the model needed 
for the nonlinear analysis of a structure through (1) the decomposition of the asym- 
metric load into symmetric and antisymmetric components and (2) the use of the 
multiple-parameter reduction technique described in the preceding section. For con- 
venience, the procedure is described herein for the case of a toroidal shell with a 
symmetric meridian (mirror symmetry in a plane normal to the axis of revolution of the 
toroid) subjected to asymmetric ring load (which is independent of the circumferential 
coordinate 0) . (See fig. 2.) The meridian is assumed to be in the form of an 
incomplete ellipse. 

The key elements of the procedure are as follows: 

1. The given asymmetric ring load is decomposed into symmetric and antisymmetric 
systems of ring loads; that is. 


P. = 
i 


( 1 ) .( 1 ) 


+ q 


( 2 ) 


p (2) 


(29) 


where q, q^, and q^ are normalizing load parameters; P^ are the normalized 
components of the given asymmetric load; and pf 1 ^ and pf 2 ^ are the normalized 
components of the symmetric and antisymmetric loads. (See fig. 2.) 

2. The basis vectors (path derivatives with respect to q^ and q^) are 

evaluated at q^=q^=o. Only hat fa the a hell moxidAjxn -Li ti&ed tn gcncAating 
tkc&C V2.ct.0-U, and the following symmetry (and antisymmetry) conditions are used for 
the generalized displacements: 
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(30) 


where Ai is in this case equal to the symmetric load parameter and Ao is 

/ o ) ^ 

equal to the antisymmetric load parameter q v 1 . Note that the symmetry conditions 
at the centerline are independent of ffl and that there are two sets of conditions 
corresponding to even and odd values of n. 


3. The reduced equations (for the full toroidal shell) are generated and the 
nonlinear response is traced by marching with the reduced equations in the two- 

dimensional load space spanned by the load parameters q^ and q^ . Each pair 

of values (q^ and q( 2 *) is selected to correspond to a value of the given load 
parameter q. (See eq. (29).) 
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The following four comments with regard to the computational procedure are in 
order : 

1. Although the individual basis vectors exhibit symmetry (or antisymmetry), 

their linear combination which approximates the solution for the original 
asymmetric load is generally asymmetric. 

2. The reduced equations can be formed by using only half the shell meridian. 

However, the solution of these equations approximates the asymmetric 
response over the full meridian. 

3. In case the given load is antisymmetric, the procedure can be simplified since 

only one X is needed. Obviously, for symmetric loads, the response is 
symmetric - 

4. Updated sets of basis vectors may be formed by using the procedure, but it is 

essential that the shell configuration used in generating these vectors be 
symmetric (e.g., q^ = 0). 

NUMERICAL RESULTS 

Accuracy of Mixed and Displacement Models 

To test and evaluate the performance of the mixed models with discontinuous 
stress resultants developed herein, several large-rotation, nonlinear problems of 
shells of revolution have been solved with these models. In all the problems con- 
sidered, the material is assumed to be linearly elastic. Comparisons are made with 
converged mixed and displacement model solutions and with previously published solu- 
tions (whenever available) to assess the accuracy and convergence of different models. 
The results of the following three problem sets are presented herein: (1) isotropic 

circular toroid subjected to uniform internal and external pressure loads; (2) iso- 
tropic spherical cap subjected to a concentrated load at the apex; and (3) laminated 
anisotropic elliptical toroid subjected to combined internal pressure and ring loads. 

Because of the axial symmetry of the load and the shell response, only one 
meridian is considered and analyzed by using the following: (1) the mixed models 

developed herein (MDm-m and MDm-n) ; (2) the standard displacement model based on 

normal integration (DEm) ; and (3) the displacement models based on reduced integration 
(DRm-n) . The characteristics of the different finite-element models are summarized in 
table 1. For all the problems considered the mixed models MD2-1 and MD3-2 are equiva- 
lent (in the sense described in the section entitled "Equivalence Between Mixed Models 
and Displacement Models") to the corresponding reduced-integration displacement models 
DR2-1 and DR3-2. The toroidal shells considered herein have symmetric meridians, and 
therefore only half the meridians of these shells are analyzed. In all three prob- 
lems, an arc-length control method is used (see refs. 9 and 11) to limit the load 
step size and to circumvent the difficulties associated with the singularity of the 
stiffness matrix near critical points. 

A solution based on the multisegment integration technique for the first problem 
is presented in reference 13. Finite-element solutions for the same problem are given 
in references 14 and 15. Analytic solutions for the spherical cap problem (for 
slightly different data) are given in reference 16. Finite-element solutions for the 
spherical cap problem are given in reference 17. 
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Isotropic circular toroid subjected to uniform internal and external pressure 
loads . - The first problem considered is that of an isotropic circular toroid subjected 
to uniform internal and external pressure loads p Q . The material and geometric char- 
acteristics of the shell and the responses of the shell to uniform internal and 
external pressures are shown in figure 3. Also shown are the actual deformed merid- 
ians of the shell corresponding to three different values of the external pressure p Q 
(designated by 1, 2, and 3). The response of the shell to external pressure exhibits 
limit-point behavior. After the first limit point, the displacements and strain 
energies of the shell are considerably larger for external pressure than for the same 
numerical values of internal pressure. Previous solutions to this problem were 
limited to the case of external pressure with magnitude of |p 0 / E | - 1.1 x 10“^ (no 
limit point) . 

For all the loads considered, the transverse shear strain energy is negligible. 
The extensional energy is the dominant energy (more than 50 percent of the total 
energy) for internal pressure load and for external pressure up to the limit point. 

The bending energy builds up rapidly after the limit point and becomes the dominant 
energy at |w c /h| ^ 2.0. It reaches 86 percent of the total energy at |w c /h| = 12.0. 
For the loads considered, the predictions of the moderate-rotation Sander s-Budiansky 
theory (refs. 18 and 19) are almost identical to those of Reissner's large-rotation 
theory used in the present study. (See fig. 3.) Also, the response to displacement- 
dependent follower load (live load) is nearly identical to that of a constant- 
directional normal load. 

The accuracy and convergence of the solutions obtained with the different models 
listed in table 1 are indicated in figures 4 and 5. The standard of comparison is 
taken to be the solution obtained by using 24 MD4-3 elements. Figure 4 shows the 
accuracy of the displacement w c /h and the total strain energy of the shell obtained 
with the different displacement and mixed models. Figure 5 shows the distribution of 
the displacements and stress resultants obtained with the different mixed models along 
the shell meridian at p Q /E = “1-1 x 10“^. As expected, the solutions obtained by 
using 24 DE2 elements are far removed from the converged solution. In the load range 
considered, no limit point is predicted with this model. The accuracy of the solu- 
tions improves by using 12 DE3 elements. Further improvement is obtained by using 
eight DE4 elements (all having the same total number of degrees of freedom). On the 
other hand, the solutions obtained with the present mixed models are considerably 
more accurate than those obtained with the corresponding displacement models. For the 
same total number of degrees of freedom, the predictions of the MD4-3 model are more 
accurate than those obtained with the MD3-2 and MD2-1 models. This is particularly 
true in the region along the meridian where sharp variations occur in the various dis- 
placements and stress resultants. (See fig. 5.) 

Isotropic spherical cap subjected to concentrated load at the apex .- The second 
problem considered is that of the nonlinear response of an isotropic spherical cap 
subjected to a concentrated load at the apex. The concentrated load retains its 
vertical direction throughout the deformation process. The shell properties are 
given in figure 6. The cap is simply supported at its boundary. The response of the 
shell exhibits both snap-through as well as snap-back phenomena. The transverse 
shear strain energy is negligible throughout the load range considered. The nonlinear 
responses predicted with the present large-rotation theory and the moderate-rotation 
Sanders-Budiansky theory are shown in figure 6. Also, the deformed configurations of 
the shell for five different values of concentrated load, as predicted with the large- 
rotation theory, are shown in figure 6. The five loads correspond to points 1 to 5 in 
the lower left-hand part of the figure. As can be seen in figure 6, the agreement 
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between the predictions of the two theories is limited to the initial phase of the 
response (with maximum displacements of the order of 4. Oh and maximum rotation of 
the order of 30 ) , after which the moderate-rotation theory overestimates both the 
displacements and the total strain energy of the shell. 

In the load range considered, the bending energy is the dominant energy. At 
q = PR/Eh = -3.56, the bending energy is 87.3 percent of the total strain energy. 

The extensional and transverse shear strain energies are only 9.7 percent and 3.0 per- 
cent of the total strain energy. 

The accuracy of the maximum displacement w c and of the total strain energy u 
obtained by using 12 MD2-1 elements, 6 MD3-2 elements, and 4 MD4-3 elements (all hav- 
ing the same total number of degrees of freedom) is indicated in figure 7. As can be 
seen in figure 7, only a slight improvement in accuracy (at higher loads) is obtained 
by using the higher order mixed models for the same total number of degrees of 
freedom . 

Laminated anisotropic elliptical toroid subjected to combined internal pressure 
and ring loads .- As a final example, consider the nonlinear response of the 10-layered, 
anisotropic elliptical toroid shell shown in figure 8. The shell has clamped edges 
(u = w = <j> = 0) and is subjected to combined internal pressure p Q and a ring 
load P. The nonlinear responses predicted with the present large-rotation theory 
and with the moderate-rotation Sanders-Budiansky theory are shown in figure 9. To 
simplify the analysis, the anisotropic bending-extensional coupling coefficients of 
the shell laminate are neglected (i.e., the shell is treated as orthotropic). Two 
cases of internal pressure loads are considered. In the first case, the internal 
pressure has an intensity p Q per unit area of the undeformed middle surface and 
is assumed to remain normal to the undeformed middle surface throughout the deforma- 
tion process. The second case is that of a live load of intensity p per unit 
area of the deforming middle surface, and the load remains normal to the deformed 
middle surface. Also, the deformed configurations of the shell corresponding to 
P q /E l = 0.66 x 10 ^ and to three different values of the ring load P as predicted 
with the large-rotation theory are shown in figure 10. The three ring load values 
correspond to points 1, 2, and 3 in the left figure. Note that for shear deformation 
shell theory, the clamped boundary conditions used do not imply that dw/ds =0. As 
can be seen in figure 9, the predictions of the moderate-rotation theory are con- 
siderably in error, particularly for higher loads. Also, the response of the shell 
to live loads is substantially different from the response to constant, directional 
loads . 

The transverse shear strain energy is significant. For example, at 

q = 0.66 x 10 and q^ ^ = -11.6, the transverse shear strain energy is 43 per- 
cent of the total strain energy. (The shear correction factor in the constitutive 
relations is taken to be 5/6.) The respective extensional and bending energies at 
the same values of q^ and q^) are 21 percent and 6 percent of the total strain 
energy . 

The accuracy of the maximum radial displacement u c . and of the total strain 
energy U obtained by using six MD2-1 elements, three MD3-2 elements, and two MD4-3 
elements in half the meridian (§ = 0 to £ = 1.0) is indicated in figure 11. In this 
case, for the same total number of degrees of freedom, the accuracy of the mixed 
models is insensitive to the degree of interpolation functions used. 
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Accuracy of Degrees-of-Freedom Reduction Technique 

In order to assess the accuracy and effectiveness of the reduction method when 
applied in conjunction with the mixed models, the isotropic circular toroid considered 
in a preceding subsection is analyzed with this method. Typical results are presented 
in figure 12. The structure is modeled by using 24 MD3-2 elements in half the meridian 
(with a total of 143 nonzero displacement degrees of freedom and 288 stress-resultant 
degrees of freedom) . Nine basis vectors (path derivatives with respect to a gener- 
alized arc length in_the solution space; see ref. 9) are generated for the unloaded 
structure (q = 0, {h} = 0, and {X} = 0) and are thus obtained by solving a linear 

set of finite-element equations. (See appendix B.) The basis vectors are ortho- 
normalized with the Gram-Schmidt procedure in order to improve the conditioning of 
the reduced equations. As a measure of the error of the reduced equations, a weighted 
Euclidean norm of the residual vectors {R} and {R} is introduced, namely 

e = 1 J{R} t {R} l {R^ {R} (31) 

q (64 + 3m) v 

where 

(R) = — [F ] {H} + [$] {x} + {© (X) } - m (32) 

{R} = [5 fc ] {H} + {TO (H»X) } - q (Q(X)} (33) 


The terms 64 and 3m are the total number of stress-resultant parameters and dis- 
placement degrees of freedom in the model. The error norm e is monitored and the 
basis vectors are updated when e exceeds the prescribed tolerance of 0.005. Nine 
new basis vectors are generated at w c /h = -3.29. The new basis vectors are used to 
advance the solution up to w c /h = -12.0. The error norm is less than the prescribed 
tolerance up to that value of w c . 

As can be seen in figure 12, the agreement between the solutions with the full 
system and the reduction technique is remarkably good except in the region just 
before the point of updating the vectors. (Updating of basis vectors is indicated 
by dotted lines between points.) Accurate solutions in that region can be obtained 
by using the updated (new) basis vectors through backtracking the solution path. At 
w c /h = -12.0 and Ua/Eh 4 = 320.0, the errors in the displacement w c and in the 
total strain energy U are 0.1 percent and 0.3 percent. 


Accuracy of Model Size Reduction Technique 

As a demonstration of the effectiveness of the procedure outlined in a previous 
section for reducing the size of the analysis model, consider the laminated aniso- 
tropic elliptical toroid shell shown in figure 8. The shell has clamped edges and is 
subjected to combined hydrostatic pressure (live load) and asymmetric ring load. 

(See fig. 13.) The hydrostatic pressure is uniform and has an intensity P 0 /Ep of 
0.66 x 10“3 p er unit area of the dd^Ohm^inQ middle surface of the shell and remains 
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normal to that, surface. The ring load is normal to the undefiomed middle. AUA^ClCt of 
the shell and is assumed to retain it* original direction throughout the analytic. 

The response of the elliptical toroid shell along with the deformed configurations at 
three different values of the ring load are shown in figure 13. The shear correction 
factors in the constitutive relations are taken to be 5/6, and the transverse shear 
strain energy becomes the dominant energy (over 50 percent of the total strain energy) 
at w c/ h = and reaches 90 percent of the total strain energy at w c /h = -5.0. 

The asymmetric ring load is decomposed into symmetric and antisymmetric compo- 
nents as shown in figure 2. The two path parameters A n and Ao are selected to be 

(1) (2) I 2 

the two load parameters q and q associated with the symmetric and antisym- 
metric ring load components. Ten basis vectors are generated at zero values of Aj_ 
and A 2 (but at uniform nonzero p Q , resulting^in a symmetric configuration) . The 
basis vectors include the nonlinear solution ({H} and {x}) evaluated at 
q = P a/E L h 3 = 0, q^ 1 ) = q^ = Q, and P Q /E L = 0.66 x 10“ 3 and all its first, 

second, and third derivatives with respect to A]_ and A 2 . These basis vectors and 
the reduced equations were generated by UAing only hal the 4 hell meridian, and the 
symmetry conditions (eqs. (30)) were applied at the centerline. 

The accuracy of the solutions obtained with the foregoing procedure for the 
range 0 - |q| 5 28.0 is indicated in figure 14. Note that the same set of basis 
vectors is used throughout the analysis, and the error norm is found to be less than 
the prescribed tolerance of 0.005. At q = -26.3, the maximum errors in the displace- 
ment w c and the total strain energy u are 0.9 percent and 2.0 percent. 

For comparison, the single-parameter reduction method is applied to this problem. 
Four, five, and six basis vectors (path derivatives with respect to q) are generated 

at q = 0 (p q /e l = 0.66 x io~ 3 ) . The generation ofi thet>e basis vectors and. the 
associated reduced equations require the use o{) the entire shell meridian. The 

accuracy of the predictions of the single-parameter reduction method with four, five, 
and six vectors is given in figure 14. As can be seen in figure 14, the solutions 
obtained with four vectors are inaccurate for |q| > 13.0. On the other hand, the 
solutions with five and six vectors are highly accurate throughout the load range 
considered. Since the major cost of the analysis is that of generating the basis 
vectors and the reduced equations , the two-parameter reduction method is computa- 
tionally more efficient than that of the single-parameter method. 


DISCUSSION OF COMPUTATIONAL STRATEGY 

The numerical results presented in this paper clearly demonstrate the accuracy 
and effectiveness of the proposed computational strategy for predicting the nonlinear, 
large-rotation, axisymmetric response of tires. In particular, the following comments 
regarding the mixed models and reduction technique are in order: 

1. In spite of the equivalence discussed in the section entitled Equivalence 
Between Mixed Models and Displacement Models," the mixed models developed herein com- 
bine the following advantages over their equivalent displacement models: 

A. The development of mixed models is simpler and more straightforward than 
that of the displacement models. This is particularly true for large- 
rotation and large-strain problems for which the functional of the mixed 
variational principle is simpler than that of the minimum potential energy 
principle . 
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B. The evaluation of the element characteristic arrays for the mixed models 


and their combinations to form the arrays K. 

, 1 D 1 

n a 'P (and their counterparts with bars and tildes) 


n (a,3) 


1 D 


and 


v (a, 3) 

r ijk 


A. (a) 

1 


B. (a) 


c (a) 

ijk' 


(see eqs . (11) to (20)) involves fewer arithmetic operations than the formu- 

lation of the corresponding arrays in the displacement models. Moreover, 
the effort required in evaluating the aforementioned arrays can be reduced 
by selecting the interpolation functions for the stress resultants to be 
OAthonOAmaZ With Aspect to the. weighting function r, thereby simplifying 

the generation of (F (See appendix B.) 


C. The elimination of the stress-resultant parameters Hj on the element 
level (condensation process) is similar to that used in hybrid models in 
which the stress-resultant field is described within the element and an 
independent displacement field is defined on the element boundaries. How- 
ever, the present mixed models differ from hybrid models in the fact that 
both stsieSA-fies ultant and dls placement fields asm described within the 
element. 


D. If the stress resultants are eliminated on the element level, the present 
elements can be viewed by a USeJl OS displacement models and can be easily 
combined with other types of displacement elements to model a shell. 

2. When reduction methods are used in conjunction with mixed models, the follow- 
ing two major advantages can be identified: 

A. The arrays appearing in the finite-element equations of the mixed models 
have a slmpleA ^0H,m than the corresponding arrays of the displacement 
models, and the evaluation of the path derivatives involves fewer arith- 
metic operations than that of the corresponding displacement models. 

b. The accuracy ofi the deduction method when used In conjunction with the 
mixed models Is expected to be hlgheA than that of the corresponding dis- 
placement method (based on using the same number of basis vectors) . This 
is particularly true for the stress resultants and is attributed to the 
fact that the path derivatives of both the stress resultants and the dis- 
placements are used as basis vectors. As a consequence of this, the basis 
vectors in the mixed method are less frequently updated than those in the 
displacement method. 

3. The procedure outlined for reducing the size of the finite-element model 

of a symmetric shell through the decomposition of the asymmetric load into symmetric 
and antisymmetric components coupled with the use of the multiple-parameter reduc- 
tion technique can lead to more dramatic savings in computational effort for 
nonaxisymmetric-response studies. Moreover, the same procedure can be applied 
when the initial discretization is done with numerical and approximation techniques 
rather than with finite elements. 


CONCLUDING REMARKS 

An effective computational strategy has been presented for the large-rotation 
nonlinear analysis of shells of revolution with application to tires. A total 
Lagrangian description of the shell deformation is used, and the analytical 
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formulation is based on a form of Reissner's large-deformation theory, including the 
effects of transverse shear deformation, laminated orthotropic material response, and 
moments acting about the normal to the middle surface. Only axisymmetric deformations 
of the shell are considered, and the load can be displacement dependent as well as 
nonconservative. The generalized stiffness coefficients (element characteristic 
arrays) are obtained by using the Hellinger-Reissner mixed variational principle. 

The polynomial interpolation (or shape) functions used for approximating the stress 
resultants are of lower degree than those used for approximating the generalized dis- 
placements. The three key components of the presented computational strategy are as 
follows: (1) use of mixed finite-element models with discontinuous stress resultants 

at the element interfaces, thereby allowing the elimination of these stress resultants 
on the element level; (2) considerable reduction in the total number of degrees of 
freedom through the use of a multiple-parameter reduction technique based on pertur- 
bation expansion; and (3) reduction in the size of the analysis model through the 
decomposition of asymmetric loads into symmetric and antisymmetric components coupled 
with the use of the multiple-parameter reduction technique. 

Numerical results have been presented for toroidal and spherical shells subjected 
to constant -directional and follower forces (live loads) . These examples demonstrate 
the high accuracy of the mixed models and show the effectiveness of using the pre- 
sented computational strategy with these models. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
August 1, 1984 
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APPENDIX A 


FUNDAMENTAL EQUATIONS OF THE NONLINEAR THEORY 
OF SHELLS OF REVOLUTION 

The fundamental equations of Reissner's large-deformation theory for shells of 
revolution are given in reference 2 and are summarized herein. Only axisymmetric 
deformations are considered, and the effects of laminated orthotropic material 
response, transverse shear deformation, and moments turning around the normal to the 
middle surface are included. A total Lagrangian description of the shell deformation 
is used, and the shell configurations at different load levels are referred to the 
initial coordinate system of the undeformed shell. 


Strain-Displacement Relationships 

The relationships between the strains and the generalized displacements of the 
middle surface of the shell are given by 


e s = cos(<f) + <(> o ) Bu - sin(<J) + <J> Q ) 8 w + cos <J) - 1 = (1 + e go ) cos y - 1 (Al) 

e 0 ■ r <“> 

k s = 3<j) (A3) 

k 0 = 7 [sin (<J> + <p Q ) - sin <j> 0 ] (A4) 

y = sin(cj) + <J> ) 3u + cos(ij) + A) 3w + sin <j> = (1 + e ) sin y (A5) 

W W O U U 

X = 1 [cos (<j> + 4> 0 ) - cos 4> 0 ] (A 6 ) 


where r _ and e n are the extensional strains in the meridional and circumferential 

• * s O v 

directions; k s and Kq are the bending strains in the meridional and circumferen- 
tial directions; y Q is the transverse shearing strain; X is the bending strain 
associated with the moment turning about the normal to the middle surface of the 
shell; e s and y are virtual extensional and shearing strains; and B = d/ds. 


Constitutive Relations 

The shell is assumed to be made of an orthotropic, linearly elastic material, and 
the relationships between the stress resultants and the strain measures of the shell 
are given by 


19 



APPENDIX A 


£ s 


a ll a l 2 | b n b 12 1 ° ° 

| | 


N 

s 

e 0 


a 22 j b 12 b 22 | ° ° 


CD 

£ 

— 


T 


— 

K s 

= 

9 11 g i2 ! 0 0 

1 


M s 

K e 


to 

to 

o 

o 


M e 

Y 


a s 0 


Q 

X 


Symm. g n 


M 

n 

— — 








(A7) 


where a a g» ^ap' (for a,P - 1 or 2), a s , and g n are shell compliance 

coefficients (inverses of shell stiffnesses) . In the present study, the stiffness 
coefficient associated with is set equal to the stiffness coefficient associated 

with Mg . Both stiffness coefficients are obtained by inverting the matrix on the 
right-hand side of equation (A7) . 
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APPENDIX B 


FORMULAS FOR COEFFICIENTS IN THE FINITE-ELEMENT 
EQUATIONS FOR INDIVIDUAL ELEMENTS 


The explicit forms of the arrays 


(a) -(a) (a) 

A I ' A I ' B Ii 


P . , and 
1 


" IJ' Ij I' I ' I ' Ij ' Ij 
are given in this appendix. For convenience, each of these 


arrays is partitioned into blocks corresponding to contributions from individual 
nodes or stress-resultant approximation functions. The expressions of the typical 
partitions (or blocks) are given in tables B1 and B2. Note that the order of the 


stress-resultant parameters in these partitions is N s , Nq, M s , Mq, Q s , and M n . 
The order of the nodal displacement parameters is u, w, and 


In tables Bl and B2 , and are respective interpolation functions for 

the stress resultants and the generalized displacements; r is the radial coordinate; 
m is the number of displacement nodes in the element; 4 is the number of parameters 
used in approximating each of the stress resultants; s^ is the element domain; 

W is the quadrature weighting coefficient (including a factor of 2tt and the 

length of the element) at the numerical quadrature point a; and 3 E d/ds. The 
range of the indices X? and j' is 1 to 4, and the range of the indices i^ 

(a) 

and j ' is 1 to m. Note that I = R . 

a 1 1 
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TABLE Bl.- EXPLICIT FORM OF TYPICAL PARTITIONS OF THE 
F, S, R, A, A, B, B, AND T ARRAYS 


Array 

Number of 
partitions 
(or blocks) 

Typical partition 

F u 

.6 x 4 

2tt f H.. r 

J s < e) * i 

a ll a i 2 j b n b 12 j ° ° 

a 22 j b 12 b 22 | ° ° 

1 - 1 _ 

*11 9 12 j 0 ° 

g 22 j 0 0 

' a s 0 

Symm. g 

^n 

ds 

s ij 

6 x m 

2t \ f U., 
J s (e) 4. 

0 0 0 

Ny o o 

0 0 r dN.' 

3 ds 

0 0 0 

0 0 0 

0 0 0 

R i 

4 

2t\ f N 

J s (e) * 

r 

0 

0 

sin <j) c 
0 

cos <t> 0 

ds 
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TABLE Bl.- Continued 
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TABLE Bl.- Concluded 


Array 

Number of 
partitions 
(or blocks) 

Typical partition 

i (a) 

ij 

£ x m 

, (a) , (a) n 

-s m <b -cos d> 0 

o o 

0 0 0 

«<«» W<“> 3N<“> r<“» 0 0 ° 

^ 3 0 0 0 

cos <j> ^ -sin <J> (a) 0 

r o r o 

0 0 0 

T (a) 

i 

m 

f o"| 

< 0 
w^> 

V J 


TABLE B2 . - EXPLICIT FORM OF TYPICAL PARTITIONS OF THE P^ AND Q i j ARRAYS 
FOR UNIFORM HYDROSTATIC PRESSURE (LIVE LOAD, q = p ) 



Number of 
partitions 
(or blocks) 

Typical partition 


m x m 

N.^ Nf . > sin <t> aW.^N.^r + N.^N... cos $ 0 

i 3 i J i 3 o 

2tt /* W. ^ 3W . r + W. ^ cos 6 0 0 ds 

J s (e) * 3 1 3 

0 0 0 

p. 

1 

m 

r sin <(> ^ 

T o 

2tt / W. ^ r ( cos cj) ) ds 

*/ s (e) 1 ] °f 

l 0 J 


24 


















APPENDIX C 


EVALUATION OF PATH DERIVATIVES 


The basis vectors are obtained by successively differentiating the governing 
finite-element equations (eqs. (24), obtained from eqs. (7) and (8)) with respect to 
the path parameters A^ and A2 « assembling the resulting equations, and solving for 
the mixed partial path derivatives. For individual finite elements, the governing 


N m+n 


equations for the path derivatives 


,m+n 


ax™ 9X2 J 


H and 


ax™ ax" 3 


X . are 


-F 


_9 H + / — X = v (m,n) 

IJ „,m ..n J ij _,m _,n j 1 


9X" ax 2 


ax 1 ax 2 


(Cl) 


,m+n 


Jl ax™ ax 2 


- h + 
n J 


(m. . 

' 13 


q (1) Q (1) 

13 


- q 


<2, a ,2,) 2 


m+n 


ax™ ax 2 3 


x. = v (m ' n) 

1 


(C2) 


where 


L Ij 


S T . + Z 
X D 


te 


) + T (d (-ta> S («) > 

D 'I ISL SLj 


COS (j) 


(a) 


- ,<«> (K 


(a) (a) > 

*1 + B IJl V 
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(a) I 


(C3) 
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(C4) 


and and are f unc tions of previously evaluated (i.e., lower order) 

path derivatives. Explicit forms for some low-order right-hand sides are as follows: 


v (1 ' 0) - v (0 ' x) 
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(C5) 
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A prescription for obtaining expressions for and is the following: 

(1) begin with -f ] . and -f ± inequations (7) and (8); (2) form their (m,n)th 

m+ \a wu-</. 

X to zero in these expressions; and 


H and 


derivatives; (3) set — 

< J ~ SA? 3 A" "i 
(4) express the results in terras of the variables A (a) , A (a) , C. (a) given by 

equations (C15) to (C20) . I I ij ' * 


In equations (C6) to (C20) , a dot (*) over a symbol refers to a derivative with 
respect to Aj., and a prime (') over a symbol refers to a derivative with respect 
to A 2 » Superscript a denotes the value of the quantity at the numerical quad- 
rature point a. The summation sign extends over all the numerical quadrature 
points. The range of uppercase indices is 1 to 64 (where 4 is the number of 
parameters used in approximating each one of the stress resultants) , and the range of 
the lowercase indices is 1 to 3m (where m is the number of displacement nodes) . 

In equations (C6) to (C12) , it is assumed that q^ = q^ = 0 . Note that the coef- 
ficients on the left-hand sides of equations (Cl) and (C2) , which must be factored, 
are the same for each of the path derivatives. Hence, the matrix of these equations 
is factored only once regardless of the number of path derivatives generated. If a 
large number of basis vectors are used, orthonormalization may be needed in order to 
improve the conditioning of the reduced system of equations - 
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If the path derivatives are generated at - 0 (and “ q^ = 0) , 

( Ci ) ( rv ) 

then H = 0 and = 0, sin <(> = <f> =0. Equations (C3) to (C20) are thus 

considerably simplified, and the computational effort in evaluating the path deriva- 
tives is reduced. For nonconservative loading, the array is asymmetric. 

However, the evaluation of the initial path derivatives (at = q^ = o) 

-involve decompo-6-vtlon ofi ci AymmeJytic. mcWvix only. 
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FORM OF THE REDUCED EQUATIONS 

The governing equations used in evaluating the undetermined coefficients \p - 
(the unknowns of the reduced equations) have the following form: ^ 
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In equations (Dl) to (D12) , the range of the indices *0, j, k, and m is 1 
to /i (the number of unknowns) ; the range of the uppercase indices is 1 to 6,6 
(where 4 is the number of parameters used in approximating each of the stress 
resultants) ; the range of the lowercase indices is 1 to 3m (where m is the number 
of displacement nodes in the element) ; and a repeated subscript in the same term 
denotes summation over its full range. Superscripts within parentheses are not 
summed . 
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SYMBOLS 






p(a) 7;(a) 
lJ ' lj 




multipliers of the trigonometric contributions to the finite- 
element equations for individual elements (see eqs. (7) 
and (8) and appendix B) 

arrays defined in equations (C15) to (C20) 


A - 0l) , A - a) , B-°J ) ,8 -°J ) , 

't 'Lj 'Lj 

j- (ot) p 

U ,K z 


a (a) n (a) ^v(a) 
i ' B ij ' C ijk' 

a' 0 ’.*.'*’.*!?’, 

1 1] 13k 

n (a,3) P (a»3) r (a,3) 

y i ' fc ij 'hjk ' 

fi(a»3) p(a#$) s(a,3) 

1 ' E ij 'hjk ' 

jj (a, 3) |(a,3) i(a,3) 

i ' lj 'ijk 


reduced-equation arrays defined in equations (D3) to (D8) 


elemental arrays (see eqs. (12) to (21) 


a 



12 ,a 22 

12 ,b 22 

12 ,g 22 


'n 


radial distance for toroidal shell (see fig. 3) 

shell compliance coefficients (inverses of shell stiffnesses) 


E Young's modulus of isotropic material 

E L ,E T elastic moduli in direction of fibers and normal to it 
e error norm 


F linear flexibility coefficients 

i j 


[?],[*] 


g lt ' g tt 

H i 

w 


global matrices defined in equations (24) 
arrays defined in equations (7) and (8) 
shear moduli in plane of fibers and normal to it 
global vectors of trigonometric contributions in <J) 
Stress-resultant parameters 

assembled vector of stress-resultant parameters 
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h 


total thickness of the shell 



linear stiffness arrays of the reduced equations (see eq, (D2) ) 

linear element coefficients defined in equation (11) 

linear stiffness coefficients for individual elements (see eq. (23)) 

generalized stiffness coefficients defined in equations (C3) and (C4) 

moment resultant turning about the normal to the middle surface of the shell 

meridional and circumferential (hoop) bending stress resultants 

external meridional concentrated moment 


m number of displacement nodes in an element 

m number of displacement nodes in the entire finite-element model 

m intensity of external meridional distributed moments 

N , N meridional and circumferential (hoop) stress resultants 

s 0 

N,N shape (or interpolation) functions used in approximating the generalized 

displacements and stress resultants 

n number of quadrature points in an element 


P ring load 

P. normalized load components 
i' i l 

p(l) p(2) normalized load coefficients of the reduced equations 

p intensity of pressure load (internal pressures have positive sign) 



intensity of distributed loads in the radial and axial directions of the 
shell 


Q transverse shear stress resultant 

s 

ofV normalized load stiffness coefficients 
13 13 

oiV oi 2) load stiffness coefficients of the reduced equations 

Q (1) ,Q (2) vectors of normalized external forces acting on the shell 

(see eqs . (24) ) 

q (1) ,q (2) load parameters 


35 



R 

R I 

m 



4 

T (a) 

i 

U 

u,w 
u, w 


radius of curvature of a spherical cap (see fig. 6) 
integrals of shape functions (see eq. (7)) 

constant vector assembled from the elemental contributions R T 
(see eqs . (24) ) 

radial coordinate of the shell 
number of basis vectors 

linear strain-displacement coefficients of the shell element (see eqs. (7) 
and (8) and appendix B) 

meridional coordinate of the shell 

number of parameters used in approximating each of the stress resultants 
within the individual elements 

number of parameters used in approximating each of the stress resultants in 
the shell 

array defined in equation (9) and appendix B 
total strain energy of the shell 

radial and axial displacement components of the middle surface of the shell 

tangential and normal displacement components of the middle surface of the 
shell 


Vi, V ± arrays defined in appendix C 

w c axial displacement at point c 

X^ nodal displacements 

{X} assembled vector of nodal displacements 

X 1 /X 2 ,x 3 Cartesian coordinate system, with x^ coinciding with the axis of revolu- 
tion of the shell 


r T ,,r . 

IaJ jsi 

[f],[r] 


matrices of basis vectors 


Y q transverse shear strain 

6 variational operator 


£ q i y virtual extensional and shearing strains 
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so 


V K e 


h' *2 

v 

V LT 

5 

n,ir 

* 

+o 


extensional strain in the meridional direction 
extensional strain in the circumferential direction 
circumferential (hoop) coordinate of the shell 
meridional and circumferential (hoop) bending strains 

bending strain associated with the moment turning about the normal to the 
shell M 

n 

path parameters 

Poisson* s ratio for isotropic material 

major Poisson* s ratio of the individual layers of a laminated shell 
nondimensional meridional coordinate 
functionals defined in equations (1) and (2) 
rotation of the shell middle surface 

= * + * 0 

angle between the axis of revolution and the normal to the shell middle 
surface 

unknowns of the reduced equations 



Ranges of 

indices : 

I,J 

1 

to 

66 

i,J 

1 

to 

66 

i, 

1 

to 

3m 

if j 

1 

to 

m 

i' , Y 

1 

to 

m 


1 

to 

/l 

S' 

1 

to 

6 
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Finite-element model notation: 


DEm 

DRm* 

MDm* 

MDm- 


displacement model based on normal integration and having m displacement 
nodes 

displacement model based on reduced integration; number of displacement 
nodes equals rri; total number of quadrature points in reduced integration 
formula equals n 


mixed model with discontinuous stress resultants at interelement boundaries 
and m displacement nodes; number of quadrature points in the Gauss- 
Legendre formula equals n; number of parameters 4 used in approximating 
each of the stress resultants within the individual elements also 
equals n 

mixed model with discontinuous stress resultants at interelement boundaries 
and m displacement nodes; number of quadrature points in the Gauss- 
Legendre formula equals m ; number of parameters 4 used in approximating 
each of the stress resultants within the individual elements also 
equals m 
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TABLE 1.- CHARACTERISTICS OF FINITE-ELEMENT MODELS USED IN 

THE PRESENT STUDY 


(a) Mixed models 


Stress- 
resultant 
appr ox ima t i on 

Displacement 

approximation 

Number of 
quadrature 
points 

Designation 

(a) 

Constant 

Linear 

1 

* 

MD2-l^ 

Linear 

Quadratic 

2 

MD3-2 

Quadratic 

Cubic 

3 

MD4-3 t 


(b) Displacement models 


Displacement 

approximation 

Number of quadrature 
points 

Designation 

(a) 

Linear 

2 

DE2 


1 

DR2-1* 

Quadratic 

3 

DE3 ** 


2 

DR3-2 

Cubic 

4 

DE4 , 


3 

DR4-3 + 


a Symbols *, ** , and + indicate equivalent finite-element 
models. 
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2 ^ 

Ua/Eh 4 

■rotatio 


essure 



X 10' 5 




(a) Two- and three-node elements. 


x i(T 5 




0 80.0 160.0 240.0 320.0 


Total strain energy, Ua/Eh 4 


(b) Four-node elements . 


Figure 4.- Accuracy of normal displacements and strain energies obtained with 
different displacement and mixed models for isotropic circular toroid 
subjected to uniform external pressure. (See fig. 3.) 
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■> 

Converged solution 


o 

MD2-1, 24 elements 

> p./E = -1.1 x 10' 5 

A 

MD3-2, 12 elements 


+ 

MD4-3, 8 elements 

J 




in-piane 

displacement 

_u_ 

h 


resultants, 

N -100 

Po h 


m 



.25 .50 

Distance, E, 


Bending 
stress 
resultants, 

M ' 20 

n h2 



.50 

Distance, E, 


Figure 5.- Accuracy of displacements and stress resultants obtained with 
different mixed models for isotropic circular toroid subjected to 
uniform external pressure. (See fig. 3.) 


fdWi\ 





Large-rotation theory 

Moderate-rotation theory 




Converged solution 

0 

MD2-1, 12 elements 

A 

MD3-2, 6 elements 

4- 

MD4-3, 4 elements 


Linear solution 


Displacement, w /h 

U 


Total strain energy, UR /Eh 


Figure 7.- Accuracy of solutions obtained with different mixed models at various 
load levels for isotropic spherical cap shown in figure 6. 





E l = 517 X 10 6 Pa (7.5 x 10 4 psi) 
E y = 8.27 x 10 6 Pa (1.2 x 10 3 psi) 
G lt = 3.10 x 10 6 Pa (450 psi) 

G n = 1.86 x 10 6 Pa (270 psi) 

V LT = ^4 

h = 1.067 x 10" 2 m (0.42 in.) 

b = 5.08 x 10" 2 m (2.0 in.) 

0 

b, = 6.223 x 10' 2 m (2.45 in.) 
b„ = 5.385 x 10" 2 m (2.12 in.) 
a = 19.558 x lo ' 2 in (7.70 in.) 


Number of layers = 10 

Fiber orientation: +45 0 /-45 0 /+45 0 /-45 0 /+45 0 /-45 0 /+45 0 /-45 0 /+45 0 /-45° 


Figure 8.- Laminated orthotropic elliptical toroid shell used in present study. 













Converged solution ^ 


o MD2-1, 6 elements 

> a (1) = — 

a MD3-2, 3 elements 

> q e l 

+ MD4-3, 2 elements 

J 










0 - 3.0 - 6.0 - 9.0 - 12.0 


Displacement, w /h 
c 



0 80.0 160.0 240.0 320.0 


Strain energy, Ua/Eh 4 

Figure 12.- Accuracy of solutions obtained with reduction method 
for isotropic circular toroid subjected to uniform pressure 
load. (See fig. 3. ) 






FiguirG 13. Deformed configurations of laminated, anisotropic elliptical toroid shell subjected to 
combined asymmetric ring load and uniform hydrostatic pressure. (See fig. 8.) 




• 28.0 




(a) Two-parameter reduction method (half model) . 



(b) Single-parameter reduction method (full model) . 


Figure 14.- Accuracy of two-parameter and single-parameter reduction methods 
for laminated, anisotropic elliptical toroid shell subjected to combined 
asymmetric ring load and uniform hydrostatic pressures. (See fig. 8.) 


53 






1. Report No. 2. Government Accession No. 

NASA TP-2343 

3. Recipient's Catalog No. 

' 4. Title and Subtitle 

MIXED MODELS AND REDUCTION TECHNIQUES FOR LARGE-ROTATION, 
NONLINEAR ANALYSIS OF SHELLS OF REVOLUTION WITH 
APPLICATION TO TIRES 

5. Report Date 

October 1984 

6. Performing Organization Code 

505-45-14-01 

7. Author(s) 

Ahmed K. Noor, Carl M. Andersen, and John A. Tanner 

8. Performing Organization Report No. 

L-15802 

10. Work Unit No. 

9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Paper 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

14. Sponsoring Agency Code 

15. Supplementary Notes 


Ahmed K. Noor: The George Washington University, Joint Institute for Advancement of 


Flight Sciences, Langley Research Center, Hampton, Virginia. 

Carl M. Andersen: The College of William and Mary, Williamsburg, Virginia; work 

performed under cooperative agreement NCCI-40. 

John A. Tanner: Langley Research Center, Hampton, Virginia. 


16. Abstract 


An effective computational strategy is presented for the large-rotation, nonlinear 
axisymmetric analysis of shells of revolution. The three key elements of the 
computational strategy are: (1) use of mixed finite-element models with discon- 

tinuous stress resultants at the element interfaces; (2) substantial reduction in 
the total number of degrees of freedom through the use of a multiple-parameter 
reduction technique; and (3) reduction in the size of the analysis model through 
the decomposition of asymmetric loads into symmetric and antisymmetric components 
coupled with the use of the multiple-parameter reduction technique. The potential 
of the proposed computational strategy is discussed. Numerical results are pre- 
sented to demonstrate the high accuracy of the mixed models developed and to show 
the potential of using the proposed computational strategy for the analysis of 
tires . 
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